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pH [ Abstract. The property of the "perfect liquid" created at RHIC is probed with 

C$ ' anisotropic flow measurements. Different initial conditions and their consequences on 

flow measurements are discussed. The collectivity is shown to be achieved fast and 
early. The thermalization is investigated with the ratio of v^jv\. Measurements from 
three sectors of soft physics (HBT, flow and strangeness) are shown to have a simple, 
linear, length scaling. Directed flow is found to be independent of system size. 
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1. Introduction: the perfect liquid 

r» , 

O . As the world's first heavy ion collider, RHIC has initiated new opportunities for studying 
\ nuclear matter under extreme conditions. After six years of successful operations, the 
discovery of the existence of a perfect liquid in ultra-relativistic heavy ion collisions 
was announcedpQ. Indications of liquid-like behavior of the matter that RHIC has 
created came in the form of large elliptic flow. Because of the pressure developed early 
in the collision, the initial spatial deformation due to geometry, which is quantified 
by eccentricity (e), is converted into the asymmetry in the momentum space, which 
is quantified by elliptic flow (1*2) [2] ■ This conversion process is directly related to the 
thermalization, equation of state, etc. The wealth of data collected and analyzed in 
many aspects, including but not limited to elliptic flow, indicates that central Au+Au 
collisions can be well described by ideal Hydrodynamics [3]. It suggests that particles in 
the medium interact with one another rather strongly, which surprised many theoretists 
who had anticipated an almost ideal, weakly interacting gas. What is more interesting 
is that, this liquid has little viscosity and acts like a perfect one [I]. This is shown in 
Fig. [H in which i>2 from data as a function of transverse momentum (jpx) is compared 
to the calculation with sound attenuation length (T s ) scaled by the time scale of the 
expansion r . The sound attenuation length is related to the shear viscosity (rj) by 
T s = |r/(e + p), where e and p are energy density and pressure, respectively. We can see 
that as expected, viscosity reduces t>2. The calculation shows that in order to explain 
the large vi observed at RHIC, one has to assume that the medium has an extremely 
small viscosity - the characteristic feature of a perfect liquid. 
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Figure 1. Elliptic flow V2 as a function of pr for different values of T s /t . The 
data points are four-particle cumulant data from the STAR Collaboration^. The 
difference between the ideal and viscous curves is linearly proportional to T s /t . This 
plot is from [4]. 
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Figure 2. Initial spatial eccentricity e at midrapidity as a function of the number 
participants for 200 GeV Au+Au collisions from various CGC models (see[7] for the 
detail description of CGC models). For comparison, the initial conditions where the 
initial parton density at midrapidity scales with the transverse density of wounded 
nucleons (full line) and of binary collisions (dotted line) are also shown. This plot is 
from [7] 



2. The initial condition 

The viscosity is so small that the initial spatial eccentricity is converted to momentum 
anisotropy with a high efficiency, and this process results in large amount of Vi as 
reported by RHIC experiments. In this explanation one assumes that the initial spatial 
eccentricity is from Glauber source[6]. Recent theoretical work (Fig. [2]) shows that a 
different initial condition like Color Glass Condensate (CGC) will give a much larger 
initial spatial eccentricity than that is from Glauber source. As a consequence of that, 
the viscosity has to be finite, as opposed to the close-to- zero viscosity in a perfect liquid, 
in order to reduce the Vi to the level that matches the data. Thus the matter that RHIC 
has created can be explained either by a perfect liquid with a Glauber source or, a viscous 
matter with a CGC source. To distinguish between the two, one has to understand the 
initial condition. However it is not easy to trace the initial condition, because with it 
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Figure 3. The directed flow v\ 
as a function of r\ for different 
Pt at b = 9 fm. Both the 
CGC model and BGK model are 
given for comparison. This plot 
is from[]. 
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the system starts, and after that the system has gone through thermalization, a possible 
QGP phase, hadronic interactions and freeze out. A lot of early information can be easily 
washed out or completely lost due to various effects at later stages. Nevertheless, both 
theoretists and experimentalists begin to realize the importance of the initial condition, 
and starts to trace its footprints. Fig. [3] shows that for high particles the v\ (solid 
lines) from CGC flips sign at r] ~ 1.2, and becomes positive for higher values of rapidity. 
That means particles are flowing in the same direction as the projective spectator. 
In the conventional factorized jet production(dashed line), the high Vi is negative 
and in the same direction as the low px bulk directed flow. It would be interesting 
for experimentalists to test this novel prediction from CGC in the future. One can also 
exam the initial condition by studying the fluctuation of elliptic flow. Both STAR[S] and 
PHOBOS [9] collaboration has measured (see FigfjJ the vi fluctuation and compared it 
to the fluctuation from initial conditions assuming Glauber sources. The V2 fluctuation 
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Figure 4. The r.m.s. width of the v-i distribution (a„2) scaled by the mean v%. Data 
are presented versus impact parameter (left by STAR) and number of participants 
(right by PHOBOS). In the left plot, together shown are eccentricity fluctuations uj (e) 
calculated from the Monte-Carlo Glauber model with standard eccentricity (crosses) 
and participant eccentricity (step- line), the latter calculation is also done by PHENEX 
(dark contour in the right plot) 

is found to be significant (~ 40% relatively), and most of it can be explained by the 
fluctuation from the Glauber model as the initial condition. It means that, again, the 
conversion process from the initial spatial eccentricity to momentum anisotropy is so 
complete that little room is left for fluctuations of other dynamic processes. 
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Figure 5. f2 as a function 
of transverse kinetic energy for 
Multi-Strange baryons. Also 
shown is the Hydro calculation 
for O and k°. The data points 
are from [IT]. 
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3. Collectivity and thermalization 

After the initial collision, particles begin to exchange momentum and the system begins 
to build up collectivity. Knowing when and how the collectivity is achieved is the first 
step towards understanding the dynamics in a hot and dense environment. This can 
be addressed by studying the v 2 of <fi and Q. Both of them are expected to have small 
hadronic cross section [10] thus are less affected by hadronic interactions. The other 
reason to choose for this purpose is because of its long lifetime - it decays outside of 
the fireball and is not formed by k + k~ coalescence, thus it picks up little information 
from a later stage. Fig. [5] shows that, although <fi and Q tends to suffer much less 
rescatterings in the hadronic stage of the collision, their Vi are found to be as high 
as other hadrons at a given px- Hence the collectivity must be developed fast at a 
pre-hardonic stage. 

Building up collectivity does not necessarily mean that the system is thermalized. 
In order for the system to be thermalized, particles in the system have to "talk" to each 
other intensively so that the information like the initial spatial anisotropy can be passed 
on to all particles. This process depends on number of collisions encountered by each 
particle. It is expected that both u 2 and are proportional to the number of collisions 
per particle, and thus the ratio of v^jv\ decreases with it [12]. In FigJHl this ratio is 
plotted against pt and compared to theoretical calculations. The Hydro calculation 
done by Borghini and Jean- Yves [12] suggests that in ideal hydrodynamics, this ratio 
decreases as a function of pt- Another version of Hydrodynamic calculation [13] shows a 
similar trend with smaller magnitude. The calculation from the AMPT[H] model shows 
a more or less flat shape. The data points are higher than theoretical calculations but 
the systematical errors are also large. It is desirable that in the future the uncertainty 
from both experiment measurement and theoretical calculation can be reduced, so that 
the degree of thermalization can be tested. 

4. Scaling of soft physics 

The number of collisions encountered by each particle on its way out not only plays an 
important role in thermalization, but also leads to a simple, but interesting scaling of 
soft physics. Fig. [7] shows that for different collision energies and over a wide range of 
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Figure 6. v^jv\ as a function 
of pt- «4 is measured by the 
three-particle cumulant method, 
and i>2 is measured by the four- 
particle cumulant method. Also 
shown are model calculations. 
This figure is from |15j 
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Figure 7. Femtoscopic radii 
dependence on the number of 
charged particles. Left panel: 
results from STAR experiment 
only for < k T >~0.20 GeV/c; 
right panel: STAR results 
combined with data from 
PHENIX,CERES and E802 
experiments, mean value of fcr 
is given on the plot. This plot is 
from 1161. 



collision systems, the HBT radii show a nice linearity if plotted against dN/dy 1 ^ 3 , which 
is proportional to the source's length, and in turn relates to the number of interactions 
for a particle on its way out. R out is an exception because it includes both space and 
time information thus the simple scaling with length is not expected. A similar dN/dy 1 ^ 3 
scaling is also observed[T7] in the strangeness yield relative to pp. Fig. [8] shows a good 
linearity if the relative yield of £1 and H are plotted as a function of dN/dy 1 ^ 3 . Also 
shown in the figure is the theoretical calculation of the enhancement with the correlation 
volume V = (N part /2) a V , where V Q = 4/3. 7iR 3 and R is the radius of the proton. The 
curve which fits the shape of the data the best is for the case of a = 1/3, which indicates 
that length plays an important role in strangeness production. Such linearity can be 
seen in flow measurements as well. In Fig. [91 the is scaled by the initial eccentricity 
and plotted as a function of paticle's density 1/SdN/dY, which is also proportional to 
the length of the system because dN/dY is proportional to the volume and S is the 
overlap area. Over a broad range of collision energies and system sizes, we observe a 
good linear relationship between Vz/e and 1/SdN/dY. This linear relation disapears 
if the same quantity plotted against N part (Fig. [TO]) , which is directly related to the 
volume. 



Collective Dynamics at RHIC 



6 




1 2 3 4 5 6 7 8 



1 2 3 4 5 6 7 8 



Figure 8. Strangeness en- 
hancement as a function of 
Pt- Also shown are three the- 
ory curves which represent the 
evolution with collision partic- 
ipants (Np ar t) of the expected 
enhancement factors. The cor- 
relation volume for strangeness 
enhancement is calculated as 
V = (N part /2) a V , where V a = 
4/3. 7ri? 3 and R is the radius of 
the proton. The three curves cor- 
respond to values of a of 1 (short 
dashed line), 2/3 (long dashed 
line) and 1/2 (solid line) respec- 
tively. The figure is re-plotted 
based on Fig. 2 in [T?] . 
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Figure 9. v 2 /e vs. 1/SdN/dy 
measured by the STAR Collabo- 
ration. Where S is the overlap 
area and dN/dy is the yield for 
charged particles at midrapidity. 
This plot is from [19]. 



Figure 10. v 2 vs. N part , for 
Cu+Cu and Au+Au collisions 
at y/s^ = 62.4 and 200 GeV. 
l-er statistical error bars are 
shown. This plot is from [20] 



The simple linear scaling from three important sectors of soft physics (HBT, 
strangeness, flow) suggests that the number of collisions encountered by each particle 
plays an important role in soft physics. One may venture [18] to predict t>2, HBT radii 
and the relative strangeness yield based on this simple scaling, without knowing anything 
about the collision (energy, system size etc.). 



Collective Dynamics at RHIC 



7 



5. Directed flow 



Directed flow (vi) describes the "bounce-off" motion of particles away from midrapidity. 
As an important tool to probe the system at forward rapidity, it complements our 
understanding of the dynamics at midrapidity. Directed flow from different energies 
at SPS has been studied in [2Tj, however its system size dependence has not been 
well explored. V\ for Au+Au collisions at both ^/snn = 62.4 and 200 GeV have been 
measured [22j. the Cu+Cu data that RHIC experiments collected in year 2005 at the 
same two energies gives us a good opportunity to study the system size dependence. The 
left plot in Fig. [TT] presents v\ as a function of pseudorapidity measured by the STAR 
Collaboration. Data from Cu+Cu collisions and Au+Au collisions at both energies 
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Figure 11. Left: Chargcd-hardon Vi vs. ry, for Au+Au and Cu+Cu collisions at 
y/sNN = 62.4 and 200 GeV. Right: The same data but plotted as a function of v\ vs. 
f] — Ubearn- Both plots arc from |23j . 

(•\/^nn — 62.4 and 200 GeV) are shown. The data points fall into two bands, one 
is for v /snn = 62.4 GeV and the other one is for ^/snn = 200 GeV. From Au+Au 
collision to Cu+Cu collision the system size is reduced by 1/3, however the V\ does 
not change. This is true even for the region near midrapidity, where t> 2 for Cu+Cu 
collisions is considerably lower that that for Au+Au collisions [TJ5]. Unlike v 2 /e which 
scales with system length, v\ is found to be independent of system size. Instead, it 
scales with the incident energy. A possible explanation to the different scalings for 
t^/e and v\ might comes from the way in which they are developed : To produce 
intensive momentum exchanges among particles are needed (and remember number of 
momentum exchanges is related to the length), while to produce V\, one in principle 
needs only different rapidity losses, which has a connection to the incident energy, for 
particles having different distances away from the central point of the collision. 

One may also test the limiting fragmentation hypothesis [24]. which has successfully 
described particle's yield and flow at forward rapidity, with different system sizes. The 
right plot in FigJH] re-plotted the same v\ results as a function of t] — ybeam- We can 
see that within three units from beam rapidity, most data points fall into a universal 
curve. This extends the validity of limiting fragmentation to different collision system 
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sizes. There are evidences |15] show that limiting fragmentation also works for higher 
harmonics like U4. 

6. Summary 

In summary, rich results from RHIC support a Hydro dynamic expansion of a thermalized 
fluid, in which the collectivity is achieved fast and at the very early time. Understanding 
the initial condition plays a key role in understanding what happens thereafter. Studying 
elliptic flow fluctuation, as well as directed flow for high p? particles, may help us 
constraint the initial condition. A few key observables from soft physics are found scaling 
with system length, which is directly related to the average number of interactions for a 
particle on its way out. Directed flow is found to depend on the incident energy but not 
on the system size. Limiting fragmentation holds for different collision energies, systems 
and flow harmonics. 
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